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Abstract 

Recent work in the literature has shown that the one-loop long distance quantum corrections to 
the Newtonian potential imply tiny but observable effects in the restricted three-body problem of 
celestial mechanics, i.e., at the Lagrangian libration points of stable equilibrium the planetoid is 
not exactly at equal distance from the two bodies of large mass, but the Newtonian values of its 
coordinates are changed by a few millimeters in the Earth-Moon system. First, we assess such a 
theoretical calculation by exploiting the full theory of the quintic equation, i.e., its reduction to 
Bring-Jerrard form and the resulting expression of roots in terms of generalized hypergeometric 
functions. By performing the numerical analysis of the exact formulas for the roots, we confirm 
and slightly improve the theoretical evaluation of quantum corrected coordinates of Lagrangian 
libration points of stable equilibrium. Second, we prove in detail that also for collinear Lagrangian 
points the quantum corrections are of the same order of magnitude in the Earth-Moon system. 
Third, we discuss the prospects to measure, with the help of laser ranging, the above departure 
from the equilateral triangle picture, which is a challenging task. On the other hand, a modern 
version of the planetoid is the solar sail, and much progress has been made, in recent years, on 
the displaced periodic orbits of solar sails at all libration points, both stable and unstable. The 
present paper investigates therefore, eventually, a restricted three-body problem involving Earth, 
Moon and a solar sail. By taking into account the one-loop quantum corrections to the Newtonian 
potential, displaced periodic orbits of the solar sail at libration points are again found to exist. 
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I. INTRODUCTION 


From the point of view of modern theoretical physics, the logical need for a quantum 
theory of gravity is suggested by the Einstein equations themselves, which tell us that gravity 
couples to T mi/ , the energy-momentum tensor of matter, in a diffeomorphism-invariant way, 
by virtue of the tensor equations QQ 


R 




1 o_ 87rG T 
v IX ' 4 ± W 


( 1 . 1 ) 


When Einstein arrived at these equations, although he had already understood that the 
classical Maxwell theory of electromagnetic phenomena is not valid in all circumstances, the 
only known forms of v were classical, e.g., the energy-momentum tensor of a relativistic 
fluid, or even just the case of vacuum Einstein equations, for which u vanishes. In due 
course, it was realized that matter fields are quantum fields in the first place (e.g., a massive 
Dirac field, or spinor electrodynamics). The quantum fields are operator-valued distribu- 


s. 


tions [3J|, for which a regularization and renormalization procedure is necessary and even 
fruitful. However, the mere replacement of by its regularized and renormalized form 
(7 ), u ) on the right-hand side of Eq. (1.1) leads to a hybrid scheme, because the classical 
Einstein tensor R^ u — ^g^R is affected by the coupling to (T^ v ). The question then arises 
whether the appropriate, full quantum theory of gravity should have field-theoretical nature 


or should involve, instead, other structures (e.g., strings |4j or loops 


and at least 16 respectable approaches [j] have been developed so far in the literature. To 
make such theories truly physical, their predictions should be checked against observations. 
For example, applications of the covariant theory lead to detailed predictions for the cross 
sections of various scattering processes js], but such phenomena (if any) occur at energy 
scales inaccessible to observations, and also the effects of Planck-scale physics on cosmology, 
e.g., the cosmic microwave background radiation and its anisotropy spectrum jli-lS], are not 
yet easily accessible to observations, although cosmology offers possibly the best chances for 
testing quantum gravity jl^ |. 


51 or twistors 


0 ). 


Recently, inspired by the work in Refs. 


15 


22] on effective field theories of gravity, where 


it is shown that the leading (i.e., one-loop) long distance quantum corrections to the New¬ 
tonian potential are entirely ruled by the Einstein-Hilbert part of the full action functional, 


some of us 


23 


24] have assumed that such a theoretical analysis can be applied to the 


3 

















long distances and macroscopic bodies occurring in celestial mechanics 25H27I]. More pre¬ 
cisely, the Newtonian potential between two bodies of masses m A and ms receives quantum 
corrections leading to 22] 


K}M = 


GmAiriB 


k\ k 2 
1 H-1—y 


0(G 2 ), 


where 


23] 


k\ = k i 


G(m A + m B ) 


/ n \2 

fc 2 = ft 2 — = «2(tp) • 
c 3 

Equation (1.2) implies that, Ve > 0, there exists an r 0 value of r such that 

GrriArriB 


V Q (r 


k\ k 2 

1 H-1- t: 


< e, Vr > r 0 . 


( 1 . 2 ) 

(1.3) 

(1.4) 

(1.5) 


This feature will play an important role in our concluding remarks in Sec. V. 

We also stress that the dimensionless parameter K\ depends on the dimensionless parame¬ 
ter k 2 . In other words, k\ is a post-Newtonian term which only depends on classical physical 
constants, but its weight, expressed by the real number ki, is affected by the calculational 
procedure leading to the fully quantum term k 2 , where the real number n 2 weighs the Planck 
length squared. More precisely, the perturbative expansion involves only integer powers of 
Newton’s constant G: 


V Q (r) ~ 


GrriArriB 

r 


OO 

E/pMG” 

p =1 


GrriArriB 

r^j - 

r 



( 1 . 6 ) 


where, upon denoting by L A and L B the gravitational radii L A = ^3—, L B = one has 
the coefficients k n = k n (L A + L B , (Ip) 2 ). At one loop, i.e., to linear order in G, where 


(m A + m B ) 1 , hi 
fi(r) = ki - - -h k 2 


q2 q 3 y>2 


(1.7) 


we can only have the contribution i - LA + LB ' > with weight equal to the real number K \, and 
the contribution with weight equal to the real number k 2 . Although the term is 
overwhelmed by the term {IjA + l G , ^he ^ wo are inextricably intertwined because is not a 
free real parameter but depends on k 2 . both K\ and k 2 result from loop diagrams. Thus, 
the one-loop long distance quantum correction is the whole term 


f / , ^2 _ (La + Lb) , (Ip) 

ji(r)G — 1—- — K\ -h n 2 —- 


( 1 . 8 ) 
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where k-i takes a certain value because there exists a nonvanishing value of k 2 . The authors 
of Ref. [22] found the numerical values 

41 


K 1 = 3, n 2 = 


107T 


(1.9) 


The work in Ref. 23] considered the application of Eqs. (1.2)-(1.4) and (1.9) to the circular 
restricted three-body problem of celestial mechanics, in which two bodies A and B of masses 
a and /3, respectively, with a > (3, move in such a way that the orbit of B relative to A is 
circular, and hence both A and B move along circular orbits around their center of mass C 
which moves in a straight line or is at rest, while a third body, the planetoid P, of mass m 
much smaller than a and /?, is subject to their gravitational attraction, and one wants to 
evaluate the motion of the planetoid. On taking rotating axes 1 with center of mass C as 
origin, distance AB denoted by /, and angular velocity u given by 


oj = 


G(a + P) 


P 


( 1 . 10 ) 


one has, with the notation in Fig. 1, that the quantum corrected effective potential for the 


circular restricted three-body problem is given by GU, where 


23] 


U = - 
2 


1 (ct + /5) 


/ 9 9 \ k 2 

(x 2 + y 2 ) + - 1 + - + -4 


s \ s s z 


( 1 . 11 ) 


73 ' iy 1 rj™ rr* 

where r and s are the distances AP and BP , respectively, while here 

G(m + a) G(m + P) 

h = «1-5-, h = %-r-, Ki = k 3 = 3. 

c c 

The equilibrium points are found by studying the gradient of U and evaluating its zeros. 


( 1 . 12 ) 


There exist indeed five zeros of gradf/ 


23] . Three of them correspond to collincar libration 


points Li, L 2 , L% of unstable equilibrium, while the remaining two describe configurations of 


stab 


23, 


e equilibrium at the points denoted by L 4 ,L 5 . The simple but nontrivial idea in Refs. 


24] was that, even though the quantum corrections in (1.2) involve small quantities, 


the analysis of stable equilibrium (to linear order in perturbations) might lead to testable 
departures from Newtonian theory, being related to the gradient of U, and to the second 
derivatives of U evaluated at the zeros of gradP. The quantum corrected Lagrange points 
L 4 and L 5 have coordinates (x(l),y + (l)) and (x(l),y-(l)), respectively, where 

{r 2 (l)-s 2 (l) + b 2 -a 2 ) 


x(l) = 


2 (a + b) 


(1.13) 


1 To be self-consistent, some minor repetition of the text in Ref. [23[ is unavoidable. 
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FIG. 1: The figure shows the two bodies of large mass, A and B, the center of mass C, and the 
planetoid at P. 


y±(l) = ±\/r 2 (l ) — x 2 (l ) — 2 ax(l) — a 2 , 


(1.14) 


where r(Z) = s(l) = w(l) ans u(l) being the real solutions of an algebraic equation 
of fifth degree (see Sec. II). Interestingly, r{l) ^ s(l), and hence to the equilateral libration 
points of Newtonian celestial mechanics there correspond points no longer exactly at vertices 
of an equilateral triangle. For the Earth-Moon-satellite system, the work in Ref. 24] has 
found 

x Q - x c ~ 8.8 mm, \y Q \ - \y c \ ~ 5.2 mm, (1.15) 

where xq (resp. xq) is the quantum corrected (resp. classical) value of x(l) in (1.13), and 
the same for yg and yg obtainable from (1.14). Remarkably, the values in (1.15) are well 


accessible to the modern astrometric techniques 


28|. 
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On the other hand, much progress has been made along the years on modern models 
of planetoids and their displaced periodic orbits at all Lagrange points Li, L 2l L 3 , L 4 , L 5 , 
to linear order in the variational equations for Newtonian theory. In particular, a mod¬ 
ern version of planetoid is a solar sail, which is propelled by reflecting solar photons and 
therefore can transform the momentum of photons into a propulsive force. Solar sailing 
technology appears as a promising form of advanced spacecraft propulsion 29l-l33j|. which 
can enable exciting new space-science mission concepts such as solar system exploration and 
deep space observation. Although solar sailing has been considered as a practical means 
of spacecraft propulsion only relatively recently, the fundamental ideas had been already 
developed towards the end of the previous century Q. 

Solar sails can also be used for highly nonKeplerian orbits, such as closed orbits displaced 
high above the ecliptic plane 351 ]. Solar sails are especially suited for such nonKeplerian 


orbits, since they can apply a propulsive force continuously. This makes it possible to 
consider some exciting and unique trajectories. In such trajectories, a sail can be used as 
a communication satellite for high latitudes. For example, the orbital plane of the sail can 
be displaced above the orbital plane of the Earth, so that the sail can stay fixed above 
the Earth at some distance, if the orbital periods are equal. Orbits around the collinear 
points of the Earth-Moon system are also of great interest because their unique positions 


are advantageous for several important applications in space mission design 


36 


|4l|. 


Over the last few dacades, severa’ 
imations of such equilibrium points 


authors have tried to determine more accurate approx- 


142 


42| . Such (quasi-)Halo orbits were first studied in Refs. 


47j. Halo orbits near the collinear libration points in the Earth-Moon system are of 


great interest, in particular around the Li and L 2 points, because of their unique positions. 
However, a linear analysis shows that the collinear libration points L\, L 2 and L 3 are of the 
type saddle x center x center, leading to an instability in their vicinity, whereas the equilateral 
equilibrium points L 4 and L 5 are stable, in that they are of the type center x center xcenter. 
Although the libration points L 4 and L 5 are naturally stable and require a small acceleration, 
the disadvantage is the longer communication path length from the lunar pole to the sail. 
If the orbit maintains visibility from Earth, a spacecraft on it (near the L 2 point) can be 
used to provide communications between the equatorial regions of the Earth and the polar 
regions of the Moon. The establishment of a bridge for radio communications is crucial 
for forthcoming space missions, which plan to use the lunar poles. Displaced nonKeplerian 
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orbits near the Earth-Moon libration points have been investigated in Refs. 29 H33L 148], 
This brief outline shows therefore that the analysis of libration points does not belong just 
to the history of celestial mechanics, but plays a crucial role in modern investigations of 
space mission design. 

Section II studies in detail the algebraic equation of fifth degree for the evaluation of 
noncollincar libration points L 4 and L 5 , by first passing to dimensionless units and then 
exploiting the rich mathematical theory of quintic equations and their roots. Section III 
derives and solves the algebraic equation of ninth degree for the evaluation of quantum 
corrections to collinear Lagrangian points L 1} L 2 ,L 3 . Section IV outlines the prospects to 
measure the quantum corrected coordinates obtained in Sec. II with the help of laser ranging. 
Section V evaluates displaced periodic orbits at the quantum-corrected Lagrange points L 4 
and L 5 , and a detailed comparison with the results of Newtonian celestial mechanics is 
also made. Concluding remarks and open problems are presented in Sec. VI, while the 
Appendices describe relevant background material on the theory of algebraic equations. 


II. ALGEBRAIC EQUATIONS OF FIFTH DEGREE FOR w(l) AND u{l) 


In Ref. 23] it has been shown that the gradU = 0 condition at noncollinear libration 
points leads to the algebraic equations of fifth degree 


where 


w 5 + C 4 W 4 + C 3 w 3 + Co = 0, 
u 5 + C+ Co — 0 , 


_ 2 K! G(m + a) _ 1 1 
C4 - 77- > Cs - ^ ~ 


2 tp 


3 k 2 c 2 lp 
7 - _ 2 «i G(m + 0) - _ ^ 

S>4 q ^ 0/2 1 ^>3 s»35 SO SO- 

O C tp 


3k 2 lpi 3 ’ 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


Such formulas tell us that it is enough to focus on Eq. (2.1), say, where, to exploit the 
mathematical theory of quintic equations, we pass to dimensionless units by defining 


1 _ 7 
w = - = 
r Ip 


(2.5) 


where 7 is a real number to be determined. The quintic equation obeyed by 7 is therefore 

7 5 + p 4 7 4 + P37 3 + Po = 0, (2.6) 












where P 4 , P 3 , po are all dimensionless and read as 


. 2 K! Ghm + a) 

Pi = (4 Ip = 77- 


Po = Co — 


3 

P3 = C 3 = 

_ 1 f Ip 
3k 2 I £ 


C 2 lp 

1 

W 

3 

= “PS 


(2.7) 

( 2 . 8 ) 

(2.9) 


At this stage, we can exploit the results of Appendix A, by virtue of which Eq. (2.6) can be 


brought into the Bring-Jerrard 


49,!50] 


form of quintic equations 


7 5 + di7 + d 0 — 0. 


( 2 . 10 ) 


Since we are going to need the roots of the quintic up to the ninth or tenth decimal digit, the 
form ( 2 . 10 ) of the quintic will turn out to be very useful, because it leads to exact formulas 
for the roots which are then evaluated numerically, which is possibly better than solving 
numerically the quintic from the beginning. Hermite 51] proved that this equation can be 


solved in terms of elliptic 


'unctions, but we use the even more manageable formulas for the 


roots displayed in Ref. (52]. For this purpose, the crucial role is played by the number 

3125 


-d n ) 4 


a = 


256 (-d^ 5 ' 

We can further simplify Eq. (2.10) by rescaling 7 according to 


( 2 . 11 ) 


The quintic for 7 is then 


7 = XX- 


~5 1 di_ do 

7 + -^ 7 + -= = 0 . 


One can choose 7 in such a way that 


- 4=1 

X 


V \" 


X = x(ii) = (—rfi) 4 . 


( 2 . 12 ) 


(2.13) 


(2.14) 


and the corresponding a of ( 2 . 11 ) reads as 

3125 ( d 0 


o = 


= a. 


(2.15) 


256 V (X(rfi)) 6 . 

If |cr| < 1, which is the case that holds in the Earth-Moon system by virtue of the numerical 
results for do and d\ obtained from the algorithm of Appendix A, then the analysis of Ref. 
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52] shows that the five roots of the quintic (2.13), here written as q 5 — 7 — j3 
obtained from the parameter 

P = 


(In 


(x(rfi )) 5 

occurring in < 7 , and from hypergeometric functions of order 4, according to 


52] 




/ i 1 Ai B 2 

4 32 32^ 


/ Fo(d)\ 

72 


1 P 5 o2 5 o3 

4 32^ 32 n 


Fi(d) 

73 


-i f -ii 7 


F 2 (d) 

w 


1 f -I/3 2 JA 3 ) 


v F s(d)y 


75 = —/ 3 Fi(a), 

where, having defined the higher hypergeometric function 


^ 1 ? ^ 2 ? •••) ttn—2 1 ) ttn—1 
bi, &2, b n - 2 , cr 


F : a —>■ F(tr) = F 
with the coefficients evaluated according to the rules 


E c - d '’ 


s=0 


n-, I S-, ( a l, s )( a 2> s)...(d n — 2 , s)(On—1, s ) 

0 ’ s (l,s)(6i,s)...(6 n _ 3 ,s)(6 re _2,s) 

(A,/r) = A(A + 1)(A + 2)...(A + p, — 1), 


one has 


Fn(a) = F 


F 1 (a) = F 


_ J _ _3_ 7 _ 11 

20 ’ 20 ’ 20 ’ 20 

113 


4’ 


2’ 4 ’ 


cr 


12 3 4 

5’ 5’ 5’ 5 
1 3 5 


.2’ 4’ 4’ 


cr 


Fold) = F 


Ffid) = F 


A 13 17 21 

20 ’ 20 ’ 20 > 20 


353 

4’ 4’ 2’ 


(7 


JL A 11 13 

10 ’ 10 ’ 10 ’ 10 

5 3 7 

4’ 2’ 4’ 


cr 


0, are 
(2.16) 


(2.17) 


(2.18) 


(2.19) 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 


This representation of the roots 7 * is discovered by pointing out that Eqs. (2.13) and (2.14) 
suggest considering such roots as functions of a = a. By taking derivatives of Eq. (2.13) 
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with respect to a up to the fourth order, one can then prove that all 7 * are particular integrals 


of the fourth-order ordinary differential equation 


52] 


ff3 ( a _ + a2 ( A ^ a ~ + a ( A2(X ~ + ( Asa ~ + C 


A = 0, 
(2.26) 


where Ai, A 2 , A 3 , B\, B 2 , B 3l C are constants. The roots 7 j undergo a peculiar variation 
when (3 describes an arbitrary curve in its plane. The critical points turn out to be 

1024 


Pi = -iC, p 2 = C,P 3 = iC, P* = -C,C = 


3125' 


(2.27) 


The group of the linear differential equation (2.26) has in this case the property that the 
root % is changed into 75 , for all k = 1,2, 3,4, when /? describes a small closed contour 
about the critical point /?*,■ 

Eventually, the roots 7 ; of Eq. (2.10) are given by 


Ai = x{di)Ai, Vz = 1,2, 3,4, 5. 


(2.28) 


At this stage, we have to invert the cubic transformation (A9) to find the five roots of the 
original quintic equation (2.6). Since in this equation the number of sign differences between 
consecutive nonvanishing coefficients is 1, we know from Descartes’ sign rule that it has only 
one positive root. We find for its numerical value (from the definition (2.5) it is clear that 
only positive values of 7 are physically admissible) 

7 + = 4.208852239482695 • 10“ 44 . (2.29) 


This value is not affected by the planetoid mass m, since m is much smaller than a in (2.7). 
As far as the unphysical roots 2 are concerned, two of them are real and negative, i.e. 

7_ = -4.20 • 10 32 or - 6.07 • 10~ 34 , (2.30) 


and two of them are complex conjugate, i.e. 

7 c = -2.10 • 10 ~ 44 ± i3.36 • 10“ 44 . 


(2.31) 


Similarly, by repeating the whole analysis for 


1 T 



2 We do not need many decimal digits for unpliysical roots. 


(2.32) 
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we find, by virtue of (2.3) and (2.4), only one positive root 

T+ = 4.208852239579132 • 10" 44 , (2.33) 

which is not affected by the planetoid mass m, since m is much smaller than /3 in (2.4), 
whereas, among the unphysical roots, two are real and negative: 

T_ = -5.17 • 10 30 or - 4.93 • 10" 32 , (2.34) 

while the remaining two are complex conjugate and read as 

T c = -2.10... • 10“ 44 ± 13.36... • 10~ 44 , (2.35) 

where the ellipsis denotes that a very tiny difference occurs in the decimal digits with respect 
to the result in (2.31), beginning at the eleventh decimal digit for the real part and at the 
tenth decimal digit for the imaginary part. 

At this stage, we can exploit Eqs. (1.13) and (1.14) to evaluate the coordinates of quantum 
corrected Lagrange points of the Earth-Moon system, finding that 

x Q = 1.8752814881103817 • 10 8 m, y Q = ±3.329001652107382 • 10 8 m. (2.36) 

The Newtonian values of such coordinates are instead 


x c = 1.8752814880224872 • 10 8 m, y c = ±3.329001652147382 • 10 8 m. 


(2.37) 


nterestingly, our detailed analysis confirms therefore the orders of magnitude found in Refs. 


23 


24j, because we obtain (cf. Eq. (1.15)) 


XQ-x c ~ 8.7894 mm, |?/q| — \yc\ ~ —4 mm. 


(2.38) 


More precisely, our refined analysis confirms to a large extent the theoretical value of xq—xc , 
whereas the sign of \yq\ — \yc\ gets reversed with respect to Eq. (1.15), and its magnitude 
gets reduced by 20 per cent. 


III. COLLINEAR LIBRATION POINTS 

From the theoretical point of view it is equally important to work out how the La- 
grangian points of unstable equilibrium, usually denoted by Li,L 2 ,L 3 , get affected by the 
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one-loop long-distance quantum corrections to the Newtonian potential. On the side of the 
applications, their importance is further strengthened, since satellites (e.g., the Wilkinson 
Microwave Anisotropy Probe) have been sent so far to the points Li,L 2 and L 3 of some 
approximate three-body configurations in the solar system. 

We beging by recalling from Ref. 23] that the gradient of the effective potential in the 
restricted three-body problem has components 


dU 

dx 


(a + 8 ) ax ( k\ k 2 \ dx f k 3 k 2 

y f -x - 5- f 1 + 2— + 3—^ ) — —5- f 1 + 2— + 3^- 


l 3 

aftl 


dU_ 

dy 


= V 


(■a + d) W 

a + d) 

~T 3 


Ri + 2 ^ + 4)-4 b+ 2 h +3 b 


S ) '^'*3 \ ryt fj- 2 


gf 1 + 2 h + 3^-4(i + 2 h + 3 h 


(3.1) 


(3.2) 


When the libration points are collinear, the coordinate y vanishes, which ensures the van¬ 
ishing of ^ as well. On the other hand, from the geometry of the problem, as shown in Fig. 


1, one has 


2 2 2 o 2 

y = r — x — Zax — a . 


The vanishing of y implies therefore that x obeys the algebraic equation 


(3.3) 


x 2 + 2 ax + a 2 — r 2 = 0, 


(3.4) 


which is solved by the two roots 


x = er — a = er — 


dl 


(a + d) 

Furthermore, the geometry of the problem yields also 

(r 2 — s 2 ) 1 (a — d) 


, e = ±1. 


x = 


21 2 (a + d) 

which implies, by comparison with Eq. (3.5), 


(3.5) 


(3.6) 


s 2 = (r — el) 2 ==>- s = ±(r — el), 


(3.7) 


where both signs should be considered, since (r — el) may be negative. Note now that the 
insertion of (3.5) into Eq. (3.1) yields 


dU d 


hi 


k. 


ae 


~ 3 ( 1 + 2— + 3 W ) - £r ) ~ — ( 1 + 2_ r + 3 W ) + 


k\ k 2 


(a + d) 
~T 3 


er 


dl 


(cr + d) 


= 0 . 
(3.8) 
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Moreover, we consider first the solution s — r — el in Eq. (3.7). This turns Eq. (3.8) into 
the form 


P 


eiy 


+ 


2 h(3 
r — eiy 


+ 


3k 2 /3 

r — eiy 


a 2 Aqa 3 k 2 a 


+ § = 0. (3.9) 


This form of the equation to be solved for r = AP suggests multiplying both sides by 
(r — el) A r 4 , which makes it clear that we end up by studying a nonic algebraic equation. 
Moreover, it is now convenient to adopt dimensionless units. For this purpose, we point out 
that the length parameters k\ and ks in the potential ( 1 . 11 ) are a linear combination of the 
gravitational radii L a , Lp of primaries and L m of planetoid according to the relations 


l ■ r, - 


Ga 


Lb = 




Lm. — 


Gm 


la L a L m: y 5 Lp T T m , 

k 1 Kil a: k$ ^3^/3 Kll/3 j 


(3.10) 

(3.11) 

(3.12) 


while k 2 = K 2 (lp) 2 from Eq. (1.4). Furthermore, all lengths involved are a fraction of the 
distance l among the primaries, and hence we set 

/ _ 1 _ P _ la _ 1/3 _ lp /q 1 

V = T, P = -1 Pa = -r, Pf3 = T , pp = -r. (3.13) 

l a l l l 

In light of (3.10)-(3.13), we find the following dimensionless form of the nonic resulting from 


Eq. (3.9): 

9 



ZA n r = o, 

72=0 

(3.14) 

where 

A 0 =— —3(1 + p) 1 n 2 (pp) 2 , 

(3.15) 


A\ = —2(1 + p) _i - 6 eK, 2 (p P ) 2 , 

(3.16) 


A 2 = — (1 + p) 1 1 — 8 eKip a + 18ac 2 (pp)" ; 

(3.17) 


A 3 = 4(1+ p) -1 £ - 3kip„ + 3eK 2 (p P ) 2 , 

(3.18) 

A 4 = 

— (1 + P) 1 (6 + (1 + e)p) — 2ki(4pq, + ppp)e + 3(1 + p)« 2 (pp) 2 , 

(3.19) 


^4-5 = (1 + p) 1 (1 + 4e) + (5 + 2e)p — 2«; 1 (p Q , + ppp) , 

(3.20) 


A(j = —(1 + p) 1 (1 + 4s) + (ICE + l)p , 

(3.21) 
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M = 2(1 + p) _1 (3 + 5p), 
^8 = — (1 + p) 1 (4 + 5p)e, 
Aq = 1. 


(3.22) 

(3.23) 

(3.24) 


If we take instead the root s = —(r — el) in Eq. (3.7) and insert it into Eq. (3.8), we find, 
with analogous procedure, the nonic equation 

9 


5^ B„r = o. 


where 


n=0 

B k = A k if k = 0,1, 2, 3,7, 8,9, 

B 4 = (1 + p) 1 (—6 + (1 — e)p) + 2«; 1 £(4p Q , + ppp) + 3 (p — 1 )k>2(pp ) 2 
-B5 = (1 + p) 1 (1 + 4e) + (5 — 2e)p — 2«i(p Q , + ppp) 

Bq = —(1 + p) 1 (1 + 4e) + (10e — l)p]. 


(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 


In Newtonian theory, the collinear Lagrangian points Li,L 2 ,L^ are ruled instead by a 
quintic equation, as is clear by setting ki — k -2 — k-s = 0 in Eq. (3.8) and multiplying the 
resulting equation by (r — el) 2 r 2 . By virtue of the two choices of sign in Eq. (3.7) one gets, 
if s — r — el, the quintic 

' S ^2 / C i^ n = 0 5 ( 2 + -„/,4 , (l + 3 P)„,.3 [1 + (1 + £ )p] „/,2 

n =0 


(1 + p) 


-£ 0 4 + 


2e 1 

-0 - 


(1 + p) (1 + p) 

while s = —(r — el) leads to the quintic 


(1 + p) 
= 0, 


- 0 3 


(1 + P) 


‘0 2 


(3.30) 


]T D n 0 n = 0 5 - 


(2 + 3p) + ^ 1 -/v .3 _ 


77.—0 


2£ 


(1+p) 


(1 + p) (1 + p) 

In this case, the coefficients are related by 


(1 + 3p) 
(1 + p) 

= 0. 


[1 - (l-e)p] , -2 


(1 + p) 


V 


(3.31) 


Ck = D k if P = 0,1, 3,4, 5, 


C 2 - — —(1 + p) 1 — 7^2,+, 


(3.32) 

(3.33) 
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C 2 , + = -(1 + p)-\l + 2 p) ^ d 2 - = -(1 + p)-\ 1 - 2 p). (3.34) 

In light of Eqs. (3.30) and (3.31), we find the values of r (distance from the planetoid to 
the primary) in Newtonian theory given by 


ri = 3.2637629578162163 • 10 8 m, (3.35) 

r 2 = 3.8167471682615924 • 10 8 nr, (3.36) 

r 3 = 4.4892055063051933 • 10 8 nr, (3.37) 

while, for the corresponding roots of the nonic equations (3.14) and (3.25), we find 

Ri = 3.263762957852764 • 10 8 nr, (3.38) 

R 2 = 3.8167471683504695 • 10 8 nr, (3.39) 

R 3 = 4.4892055063281494 • 10 8 nr. (3.40) 

By virtue of these values, we find 


Ri — r\ — 3.6 mm, 


(3.41) 


R 2 — v 2 = 9 mm, (3.42) 

R 3 — r 3 = 2.3 mm. (3.43) 

Interestingly, the order of magnitude of quantum corrections to the location of Li,L 2l L 3 
in the Earth-Moon system coincides with the order of magnitude of quantum corrections 
to L 4 , L 5 that we have found in Eq. (2.38). This may not have any practical consequence, 
since L\, L 2 , L 3 are points of unstable equilibrium, but the detailed analysis performed in this 
section adds evidence in favour of our evaluation of quantum corrections to all Lagrangian 
points in the Earth-Moon system being able to predict effects of order half a centimeter. 

The main perturbations of such a scheme may result from the Sun. If one then considers 
a restricted four-body problem where the Earth and Moon move in circular orbits around 
their center of mass, which in turn moves in a circular orbit about the Sun 3 , one finds that 
L 4 and L 5 are no longer points of stable equilibrium 53| . However, one can evaluate the 
impulse required to induce stability at L 4 , i.e. to force the planetoid to stay precisely at L 4 . 
Such an impulse turns out to be 2360 lb/sec/slug/yr, as shown in Ref. [53]. 


3 


The Sun’s effect on the planetoid is much larger than the Sun’s effect on the Moon. 


16 








IV. TINY DEPARTURE FROM THE EQUILATERAL TRIANGLE PICTURE: 
PROSPECTS TO MEASURE THE EFFECT WITH LASER RANGING 


The quantum gravity effect described in our paper can be studied with the tech¬ 
nique of Satellite/Lunar Laser Ranging (hereafter SLR/LLR) and a laser-ranged test mass 


equipped with Cube Corner Retro-reflectors (CCRs), to be designed ad hoc for t 


pose. SLR/LLR is performed by the International Laser Ranging Service (ILRS) [54J, which 


11s pur- 


recently celebrated the 50th anniversary of the first successful SLR measurement, which 
occurred at the Goddard Geophysical and Astronomical Observatory (GGAO) on October 
31, 1964 4 . Detecting this tiny departure from classical gravity is a challenging task, which 
requires precise positioning in space at the Lagrangian points L 4 and L 5 , in absolute terms, 
that is, with respect to an appropriately chosen coordinate reference system. One potential 
choice is the International Terrestrial Reference system (ITRS) (55|, which is established 
with several geodesy techniques, including SLR/LLR. The latter provides almost uniquely 
the metrological definition of the Earth’s center of mass (geocenter) and origin of the ITRS, 
as well as, together with Very Long Baseline Interferometry (VLBI), the absolute scale of 
length in space in Earth Orbit. Given its similarity with LLR 56), another option for the 
coordinate frame is the Solar System Barycenter (SSB). In fact the distance of L 4 and L5 
from the ground laser stations of the ILRS is very close to their distance to the Laser retrore- 
flector array (LRAs) of deployed by the Apollo and Lunokhod missions, which over the last 
45 years were used for some among the best precision tests of General Relativity (see Refs. 


57 


62]). The SSB is particularly apt for the purpose, since it is used for General Relativity 


tests carried out with LLR data analysis by means of the orbit software package Planetary 


Ephemeris Program (PEP) since the eighties 58| and until nowadays 59|, [60]. PEP has been 


developed by the Harvard-Smithsonian Center for Astrophysics (I. I. Shapiro et al, currently 
maintained by J. F. Chandler). A review of LLR data taking and analysis can be found in 
Ref. 60]. 

A laser ranging test mass (LR tm ) can be designed with a dedicated effort, by exploiting 
the experience of LLR data taking and analysis described above, and especially by taking 
advantage of existing capabilities for detailed pre-launch characterization of any kind of 


4 See http://ilrs.gsfc.nasa.gov for a description of satellite/lunar laser ranging. 
http://ilrs.gsfc.nasa.gov/ilrwl9/. 


See also 
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LRAs and/or test mass for Solar System exploration 163 h 66| . Some of the Key Performance 
Indicators (KPIs) that must be taken into account to design an appropriate LR tm for the 
signature of new physics described in this paper are as follows. 

(i) Adequate laser return signal (lidar optical cross section) from the Lagrangian points 
La, L§. 

(ii) Acceptable rejection of the unavoidable nongravitational perturbations (NGPs) at L 4 , L 5 
which any chosen test mass and/or test spacecraft will experience, whose complexity scales 
with the complexity of the structure of the test mass and/or test spacecraft itself. 

(iii) Optimization/ minimi zation of the value of the surface-to-mass ratio, S/M. This is a 
critical KPI, since all NGPs related to the sun radiation pressure and thermal effect, are 
proportional to S/M (see for example Ref. [67|). Compared to other test spacecrafts and/or 
test masses an LR tm has the advantage of the simplicity of geometrical shape (for example, 
spherical) and mechanical structure. To date, Apollo/Lunokhod are demonstrating a lifetime 
of at least 45 years. 

(iv) Time-durability of the test mass to prolonged measurements. Since LR tm are passive 
and maintenance free, this KPI favors LR tm over other types of any active test masses and/or 
spacecrafts. 

he above KPIs can be characterized at the dedicated laboratory described by Refs. 


63, 


64j (see also http://www.lnf.infn.it/esperimenti/etrusco/). From the experimental point 


of view of laser ranging investigations, arguments reported in this section for L 4 and L 5 apply 
identically to L 3 . They do not apply to L 2 since such a position is not visible from ILRS 
stations. The distance of L 4 from Earth is shorter than for L 3 , L 4 and L 5 , which would 
make the laser return signal from an LR tm in L 4 higher than from L 3 ,L 4 ,L 5 (by a purely 
geometric factor equal to the fourth power of the ratio of the distances of L 3 and L\ from 
any given ILRS station; see for example Ref. 63]). Given the relative proximity of L\ to 
the Moon, gravitational effects on an LR tm in L 4 related to the nonpointlike structure of the 
Moon (felt in L 4 ) should be evaluated to determine their influence, if any, on the conclusions 
of the previous section. This influence is expected to be negligible for an LR tm in L 3 , L 4 and 
L 5 , since they are much more distant from the Moon than L 4 . 
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V. DISPLACED PERIODIC ORBITS FOR A SOLAR SAIL IN THE EARTH- 
MOON SYSTEM 


Displaced periodic orbits describe the dynamics of the planetoid, e.g., a solar sail, in 
the neighborhood of the libration points, which have been studied in detail in the quantum- 
corrected case 


and in Newtonian theory 


29j. The appropriate tool of classical mechanics 


are the variational equations, for which we refer the reader to Refs. 24, [26j, [27]. In the 


simplest possible terms, the components x, y, z of the position vector of the sail (see Fig. 2) at 
each libration point change by the infinitesimal amount £,77, ( respectively and, by retaining 
only first-order terms in £, 77, £ in the equations of motion, one finds the following linear 
variational equations of motion for the libration points L 4 . L 5 describing stable equilibrium 


29] 


{ - 2 f]= + qp) + a ( , 

V + 2 £ — + a,, 

c = y°,c+a c . 


(5.1) 

(5.2) 

(5.3) 


where the auxiliary variables a^,a v ,a^ describe the solar sail acceleration, and 
Uxx , Uyy 1 U ^ z , U® y are the partial derivatives of the gravitational potential (1.7) evaluated 

29]. we are here using units where the sum of the 


at L 4 or L5. Note that, following Ref. 


masses of the primaries is set to 1, as well as their distance and the Newton constant. 


Following Ref. 


29|, assume now that a solution of the linearized equations of motion 


(5.1)-(5.3) is periodic of the form 


£(£) = A£ cos(u;*t) + B^ sin(u;*f), 
77 (f) = A v cos(ujA) + B v sm(uj*t), 


(5.4) 


(5.5) 


where A^,A V ,B^ and B v are parameters to be determined, and cn* = 0.923 is the angular 
rate of the Sun line in the corotating frame in a dimensionless synodic coordinate system 
29l ]. By substituting Eqs. (5.4)-(5.5) in the differential equations (5.1)-(5.3), we obtain the 


following linear system in A^, A^, B g and B v 


29(]: 


- K 2 + Ul)B ( + 2W.A, - = 0, 

— + 2— (wf + U° )A V — 0, 


(5.6) 

(5.7) 
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z 



FIG. 2: Schematic geometry of the Earth-Moon restricted three-body problem when the planetoid 
is a solar sail. 


- O* + UL)^ - U xy A V - 2uj * B r, = «0 COS 3 ip, (5.8) 

- 2 uj+At: - Ul y B(: - (ujI + Uyy)B v = —do cos 3 ip. (5.9) 

This linear system can be solved to find the coefficients A^, B^, A V1 B v , here arranged in 
the four rows of a column vector P, while b is the column vector whose four rows are the 
right-hand sides of (5.6), (5.7), (5.8) and (5.9), respectively. Let A be the 4x4 matrix 


A = 


where the 2x2 submatrices of A are 



(A, B , 
D r 


71 ! = 




B i = 

C\ = 



-U° \ 

0 )' 

0 V 

- u 0 / ’ 

,: v/ 


(5.10) 


(5.11) 


(5.12) 
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(5.13) 










FIG. 3: Time evolution of £ for L 4 in the Newtonian case. 


z>,= " 2w * 

0 -“l-K, 

With this matrix notation, the solution of our linear system (5.6)-(5.9) reads as [29] 


(5.14) 


P‘= (A' 1 )) b> Vi = 1,2,3,4. (5.15) 

The coefficients A^,A V ,B^ and B v are amplitudes that characterize the displaced periodic 
orbit. 


Last, the out-of-plane motion (Eq 
the solution of Eq. (5.3) is given by 


291 


4.3)) is decoupled from the in-plane motion, hence 


T o i-l 

' ZZ I 


£(t) = 8(t)ao cos 2 </?(sin <p)\U\ 

+ cos£(t = 0) — a 0 cos 2 <^(sin ip)\U\ 


o i-l 

ZZ I 


(5.16) 


where 9 is a step function 


6{t) — 1 if t > 0, 9{t) — 0 if t < 0. 


Thus, the required sail acceleration for a fixed distance can be given by 29] 

C(* = 0)|t£| 


do — 


(5.17) 


(5.18) 


cos 2 (/p(sin <p) 

In Newtonian theory, the findings for displaced periodic orbits are well summarized in 
Figs. 3-5 (cf. Ref. 29]). 

On considering the quantum corrections evaluated in detail in Sec. II, and setting fur¬ 
thermore the angle tp = j, while £(t = 0) = 2, w* = 0.923, a 0 = 10~ 4 , we arrive at the 
plots displayed in Figs. 6-8. The starting value of £ has been taken to be 100 km, increased 
gradually to reach 2500 km. 


21 












10000 


5000 


-5000 

-10000 



FIG. 4: Time evolution of r) for L 4 in the Newtonian case. 



FIG. 5: Periodic orbits at linear order around the Lagrangian point L 4 in Newtonian theory. 


Our calculation is of interest because it shows that even our quantum corrected potential 
allows for periodic solutions in the neighborhood of uniform circular motion. The precise 
characterization of regions of stability and instability [70] of such displaced periodic orbits 
is a fascinating problem for the years to come. 





FIG. 6 : Time evolution of £ for L 4 in the quantum-corrected model. 
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FIG. 7: Time evolution of r) for L4 in the quantum-corrected model. 



FIG. 8: Periodic orbits at linear order around the Lagrangian point L4. We have used the quan¬ 
tum corrected coordinates obtained at the end of Sec. II. The periodic orbit is elliptic as in the 
Newtonian case displayed in Fig. 5. 

VI. CONCLUDING REMARKS AND OPEN PROBLEMS 

We find it appropriate to begin our concluding remarks by stressing two conceptual 
aspects, which are as follows. 

(i) In the course of an orbit of a celestial body around another celestial body, their mutual 
separation may change by a nonnegligible amount. Thus, it would be misleading to look 
for an observational test of one-loop long-distance quantum corrections to the Newtonian 
potential by investigating the orbits, because we do not have a formula for V(r) which is 
equally good at all points. By contrast, the evaluation of stable equilibrium points (to first 
order in perturbations) provides a definite prediction, i.e., the coordinates of such a point, 
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which can be hopefully measured with the techniques outlined in Sec. III. In other words, 
coordinates of Lagrangian libration points L4 and L5 and displaced periodic orbits around 
unperturbed circular motion provide a valuable test of effective field theories of quantum 
gravity, whereas the orbits of celestial bodies are best studied within the framework of 
relativistic celestial mechanics. 


(ii) At the risk of repeating ourselves, the technique of Refs. 


15 


11 


22] provides corrections 


to the Newtonian potential, and hence the unperturbed dynamics is the Newtonian celestial 
mechanics of the Earth-Moon-satellite system, which may provide a good example of circular 
restricted three-body problem. The quantum corrected potential becomes (1.2), where k\ 
appears, on dimensional ground, purely classical, but includes a numerical coefficient, K \, 
which depends on the value taken by the coefficient n 2 that multiplies l 2 P in k 2 . 


Ki = Ki(k 2 ), k 2 = k 2 1 ] 


Thus, we do not compute corrections to relativistic celestial mechanics (cf. Ref. 68]), but, 
on the other hand, we need the advanced tools of relativistic celestial mechanics to test the 
tiny effect predicted in Eq. (2.38). We should mention at this stage the important work in 
Ref. |69|], where the authors obtain a triangular solution to the general relativistic three- 
body problem for general masses, and find that the post-Newtonian configuration for three 
finite masses is not always equilateral. When their technique is applied to the Earth-Moon 
system, we End, unlike our Eq. (2.38), a correction to the ^-coordinate of L4 of order 2.73 
mm, and a correction to the y coordinate of L4 of order —0.53 mm. The former agrees 
with our orders of magnitude, while the latter, being less than a millimeter, is very hardly 
detectable. 

Our original contribution is, first, the detailed calculation of the roots of the quintic 
equation in Sec. II (which is an original application of techniques previously developed by 
mathematicians), second, the derivation and solution of the nonic equation in Sec. Ill for 
quantum corrections to collinear Lagrangian points and, third, the application of the roots 
in Sec. II to the evaluation of displaced periodic orbits of solar sails in the Earth-Moon 
system, when the one-loop long distance quantum corrections to the Newtonian potentials 


are taken into account 


22 


24]. The use of dimensionless variables for the quintic equation 


(2.6), and the exploitment of exact formulas for its roots were of crucial importance to 


double-check the numerical predictions of Refs. 


23 


24] . Interestingly, we have found that 
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a refined analysis like ours confirms the orders of magnitude obtained in Ref. 24], whereas 
the sign of \i/q\ — \yc\ gets reversed with respect to Ref. 24|, and its expected theoretical 
value turns out to be smaller by 20 per cent. Furthermore, displaced periodic orbits have 
been evaluated in Sec. IV with the quantum corrected coordinates displayed in Eq. (2.36), 
when the condition for the existence of displaced orbits is affected by terms resulting from 
a solar sail model. We have found that, even when the quantum corrected potential (1.2) 
is adopted, the displaced periodic orbits are of elliptical shape (see Fig. 8) as in the New¬ 
tonian theory. The solar-sail model is an interesting possibility considered over the last few 
dacades, but is not necessarily better than alternative models of planetoid. For example, 
the large structure and optical nature of solar sails can create a considerable challenge. If 
the structure and mass distribution of the sail is complicated, one has to resort to suitable 
approximations. Furthermore, the characterization of regions of stability or instability |70| 


of displaced periodic orbits of solar sails is a theoretical problem whose solution might have 
far reaching consequences for designing space missions. 

Last, but not least, the laser ranging techniques outlined in Sec. IV appear as a promising 
tool for testing the predictions at the end of Sec. II. The years to come will hopefully tell 
us whether a laser ranging test mass can be designed, upon consideration of the four key 
performance indicators listed at the end of Sec. IV. At that stage, the task will remain to 
actually send a satellite at L 4 and keep it there despite the perturbations caused by the 
Sun 53], which has never been accomplished to the best of our knowledge. The resulting 
low-energy test of quantum gravity in the solar system would reward the considerable effort 
necessary to achieve this. 
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Appendix A: Bring-Jerrard form of quintic equations 


Let us start from the general quintic equation 


X 3 + CI4A 4 + G3 A 3 + d 2 X~ -f- d\X + CIq — 0. 


Denote the roots of Eq. (Al) by Aj, i — 1,5, and let 


(Al) 


Sn = S n (X k ) = ^(A* 


(A2) 


fc=i 


be the sum of the nth powers of such roots. By virtue of the Newton power-sum formula, a 
general representation of S n is 

n— 1 

e . 

-31 


Sn TlCl^ — n ^ ^ S n —j(l5—j 
3 = 1 


(A3) 


with the understanding that aj = 0 for j <0. For the lowest values of n, Eq. (A3) yields 
Si(X k ) = —04, S 2 (Xf c) = (0-4)^ — 203, Ss(Xk) = —(«4) 3 + 30304 — 3a2, 


— (04) 4 — ^03(04)^ + 40204 + 2 (03) 2 — 4 ai, 


(A4) 


S§(X k ) — —(04)° + 5 03(04) 3 — fl2( a 4) 2 — { a 3)^ a 4 A a 1®4 — a 0 A O2O3 . (A5) 

A systematic way to proceed involves two steps, i.e., first a quadratic Tschirnhaus transfor¬ 


mation 


III 


Vk = (Xk) 2 A ixX k + n (A6) 

between the roots X k of Eq. (Al) and the roots Y k of the principal quintic 


x + C2X + Cl x + cq — 0 , 


supplemented 72] by the evaluation of Si(Y k ),..., Ss(Y k ) to o 
fj,, v, c 0 , Ci, c 2 , and eventually a quartic Tschirnhaus transformation 


(AT) 

Dtain through radicals 

HI 


Z k — (Ifc) 4 + u i (Y k ) 3 + u 2 (Y k ) 2 + u%Y k + W4, 


(AS) 


between the roots Y k of Eq. (A7) and the roots Z k of the Bring-Jerrard form (2.10) of the 
quintic. This procedure is conceptually clear although rather lengthy (see Appendix B), 
and the joint effect of inverting (A8) and then (A6) to find X k = X k (Yj(Zi )) leads to 20 
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candidate roots [72j, which is not very helpful if one is interested in the numerical values 
of such roots, as indeed we are. One might instead stick to our Eq. (2.6) and then try 
to exploit the Birkeland theorem [73f, according to which the roots of any quintic can be 
re-expressed through generalized hypergeometric functions like the ones defined in our Sec. 
II. However, when there are three (or more) nonvanishing coefficients as in Eq. (2.6), the 
Birkeland theorem leads to too many (for numerical purposes) hypergeometric functions in 
the general expansion of roots. Furthermore, the set of linear partial differential equations 
0 obeyed by the roots when viewed as functions of all coefficients does not lead easily to 
their explicit form. 

At this stage, having appreciated the need for a trinomial, possibly Bring-Jerrard form 
of the quintic, and rather than feeling in despair, we point out that, since in our original 
quintic ( 2 . 6 ) two coefficients vanish, i.e. a 2 = a± = 0 , it is more convenient to use what is 
normally ruled out in the generic case 72], i.e., a cubic Tschirnhaus transformation between 
the roots X k of Eq. (Al) and the roots 17 of Eq. (2.10): 

Y k = (X k ) 3 + Ai(AA ) 2 + X 2 (X k ) + A 3 . (A9) 

By virtue of Eqs. (2.10) and (A2)-(A5), we find 

Si(y*) = S 2 (Y k ) = S 3 (Y k ) = 0 , (A 10 ) 

5 4 (17) = —46b, 5*5(17) = -5 d 0 . (All) 

O 11 assuming the cubic relation (A9), Eqs. (A10) become a nonlinear algebraic system 
leading to the numerical evaluation of Ai, A 2 , A 3 . More precisely, from S\ (17) = 0 we find 

5 A 3 + X 2 Si(X k ) + Ai5* 2 (A7) + S 3 (X k ) = 0, (A12) 


while from 5 * 2 ( 17 ) = 0 we obtain 

5(As) 2 + 2A 2 A35 *i(A7) + [(A 2) 2 + 2XiX 3 ]S2(X k ) + 2(AiA2 + X 3 )S 3 (X k ) 

+ [(AO 2 + 2A 2 ]5* 4 (X fe ) + 2A 1 5' 5 (A7) + S 6 (X k ) = 0. (A13) 

Last, from the vanishing of S 3 (Y k ) we get 

5(As) 3 + 3A 2 (A3) 2 5*i(A7) + 3(A 2 ) 2 As52(A7) + (X 2 ) 3 S 3 (X k ) 

+ 3(Ai) 2 A 3 5 4 (A7) + [3(Ai) 2 A 2 + 6AiA 3 ]5* 5 (A7) + [(A 1) 3 + 3 A 3 + 6AiA2]5*6(Affc) 

+ 3[(Ai) 2 + A 2 ]5V(A7) + 3XiSg(X k ) + Sg(X k ) = 0. (A14) 
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The system (A12)-(A14) cannot be solved by radicals because, if one expresses for example 
Ai as a linear function of A2 and A3 from Eq. (A12), and one solves the resulting quadratic 
equation for A 2 = A 2 (A 3 ) or A 3 = A 3 (A 2 ) from Eq. (A13), one discovers that Eq. (A14) is 
not a polynomial in A 3 (respectively, A 2 ). Nevertheless, for numerical purposes , the system 
(A12)-(A14) can be solved and has been solved from us in the Earth-Moon system. Last, 
from Eq. (All) we find the coefficients d± and do in the Bring-Jerrard form of the quintic, 
according to the formulas 

1 12 

dx = ~-S 4 (Y k ) = Y,buSi(X k ), (A15) 

i =0 

1 15 

d 0 = -- S 5 (Y k ) = b 0l S l (X k ), (A16) 

i =0 

We have evaluated all bu and b 0 i coefficients by applying patiently the Tschirnhaus trans¬ 
formation (A9) and the definition (A2). We find therefore six triplets of possible values for 
Ai, A 2 , A 3 (see Tables I and II), which lead always to the same values of d,\ and do (this is a 
crucial consistency check), i.e. 


di(w) = 2.78 • 10" 176 + i6.91 • 10” 186 , di(u) = 2.78 • 10” 176 - i8.51 • 10" 188 , (A17) 

do(w) = 1.66 • 10~ 220 - i5.17 • 10" 230 , d 0 (u) = -1.66 • 10" 220 + i6.36 • 10" 232 , (A18) 

where w and u are the variables defined in Eqs. (2.5) and (2.32), respectively. Eventually, 
the roots X k of our quintic (2.6) have been obtained by solving Eq. (A9) for X k = X k (Y k ), 
with the help of the solution algorithm for the cubic equation. This means that we first 
re-express (A9) in the form 


h(X k ) — (X k ) 3 + k 2 (AA,) 2 + K\X k + Ko — 0, 
where k 2 = A 3 , «q = A 2 , k 0 = A 3 — Y k . We then define the new variable 

v k = x k + ^ = x k + 

in terms of which Eq. (A19) is mapped into its canonical form 


(14) 3 +pV k + q = 0, p = h' 


As shown in Ref. 52], if the discriminant 
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4 pi 



(A19) 


(A20) 


(A21) 


(A22) 



is such that | <5| < 1, or if 5 = 1, the three roots of Eq. (A21) can be expressed through the 
Gauss hypergeometric function in the form 






(-1 ) 3i F 


(_1 1 1 \ 1 
V 6’6’2' ) + 3 



1 2 3 

3’ s’ 2’ d 


1 ,2), (A23) 


(V k ) 3 



12 3 
3’ 3’ 2’ ° 


(A24) 


As is clear from (A20), (A23), (A24), our method yields eventually 15 candidate roots, and 
by insertion into the original quintic (2.6) we have found the 5 effective roots at the end of 
Sec. II. 


TABLE I: The six triplets of values of A 3 , A 2 and Ai for the 1/r-equation (2.1). 


nth triplet 

A 3 

A 2 

Ai 

11=1 

-4.98 x 10 45 - 3.10 x 10 " 55 i 

0.26 + 1.59 x 10 " 11 i 

4.21 x 10 32 + 3.78 x 10 " 44 i 

n =2 

-4.98 x 10 “ 45 + 3.10 x IO - 55 i 

0.26 - 1.59 x 10 " 11 i 

4.21 x 10 32 - 3.78 x 10 " 44 i 

n=3 

2.49 x 10 -45 + 4.32 x 10 " 45 i 

0.26 + 2.57 x 10 -11 i 

4.21 x 10 32 + 6.11 x 10 -44 i 

n=4 

2.49 x IO - 45 - 4.32 x 10 " 45 i 

0.26 - 2.57 x 10 " 11 i 

4.21 x 10 32 - 6.11 x 10 " 44 i 

n=5 

2.49 x 10 “ 45 + 4.32 x 10 " 45 i 

0.26 + 9.82 x 10 " 12 i 

4.21 x 10 32 + 2.34 x 10 " 44 i 

11=6 

2.49 x 10 " 45 - 4.32 x 10 " 45 i 

0.26 - 9.82 x 10 " 12 i 

4.21 x 10 32 - 2.34 x 10 " 44 i 


TABLE II: The six triplets of values of A 3 , A 2 and Ai for the 1/s-equation (2.2). 


nth triplet 

A 3 

A 2 

Ai 

n=l 

-4.98 x 10" 45 - 3.82 x 10“ 57 i 

0.26 + 1.96 x 10" 13 i 

5.17 x 10 30 + 3.78 x 10 -44 i 

n=2 

-4.98 x 10" 45 + 3.82 x 10“ 57 i 

0.26 - 1.96 x 10-13 i 

5.17 x 10 30 - 3.78 x 10" 44 i 

n=3 

2.49 x 10" 45 + 4.32 x 10" 45 i 

0.26 + 3.16 x 10" 13 i 

5.17 x 10 30 + 6.11 x 10 -44 i 

n=4 

2.49 x 10" 45 - 4.32 x 10" 45 i 

0.26 - 3.16 x 10" 13 i 

5.17 x 10 30 - 6.11 x 10" 44 i 

n=5 

2.49 x 10" 45 + 4.32 x 10" 45 i 

0.26 + 1.21 x 10" 13 i 

5.17 x 10 30 + 2.34 x 10 -44 i 

11=6 

2.49 x 10 " 45 - 4.32 x 10 " 45 i 

0.26 - 1.21 x 10" 13 i 

5.17 x 10 30 - 2.34 x 10 " 44 i 


Yet another valuable solution algorithm is available, i.e. the method in Ref. 0 which 
expresses the roots of the quintic (Al) through two infinite series, i.e., the Jacobi norne and 
the theta series, for which fast convergence is obtained, but the need to evaluate the roots 
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with a large number of decimal digits makes it problematic, as far as we can see, to deal 


with such series. Further valuable work on the quintic can be found in Ref. 


76]. 


Appendix B: Alternative route to the quintic 

We here find it useful to give details about the main alternative to the procedure used in 
Appendix A. For that purpose, as we said, one assumes that the roots X k of Eq. (Al) are 
related to the roots Y k of the principal quintic (A7) by the quadratic transformation (A 6 ). 
The power sums for the principal quintic form are indeed 

Si(Y k ) = S 2 (Y k ) = 0, S 3 (Y k ) = —3c 2 , S 4 (W) = -4 Cl , S 5 {Y k ) = -5c 0 . (Bl) 

On the other hand, we can evaluate Si(Y k ) and S- 2 (Y k ) by using the quadratic transformation 
(A 6 ) and exploiting the identities 

Si{Y k ) = S 2 (X k ) + pS^X,) + hu, (B2) 

S 2 (Y k ) = S 4 (X k ) + 2 /iS 3 {X k ) + (p 2 + 2 u)S 2 (X k ) + 2 fivS^Xj + 5u 2 , (B3) 

obtaining therefore the following equations for p and v. 

pa 4 — 5u + 2a 3 — (a 4 ) 2 = 0, (B4) 

/i 2 a 3 — lOz / 2 + p(3 a 2 — a 3 a 4 ) + 2 a 4 — 2a 2 a 4 + (a 3 ) 2 = 0, (B5) 

where, in the course of arriving at Eq. (B5), we have re-expressed repeatedly (a 4 ) 2 from 
Eq. (B4). This system is quadratic with respect to p and u, and hence leads to two sets of 
coefficients. For the case studied in Eq. (2.6), they reduce to (here a 3 = p 3 , a 4 = p 4 ) 

a 4 [13a 3 — 4(a 4 ) 2 ] ± ^/60(a 3 ) 3 — 15(a 3 a 4 ) 2 


F± = 


2[5a 3 - 2 (a 4 ) 2 ] 

U-j- 2 1 . . o 

Rt = —04 + -o 3 - -(a 4 ) . 


(B 6 ) 


(B7) 


There is complete freedom to choose either of these. After finding p and v in such a way, 
one can use the Eqs. (Bl) to obtain c 0 ,Ci,c 2 . One finds explicitly, in general, 

c 0 = -u 5 - pv^{X k ) - ( 2 //V + v 4 )S 2 (X k ) - ( 2 //V + 4pu 3 ) S 3 (X k ) 

— (/z 4 ^ + Qp 2 v~ + 2z^ 3 ) S 4 (X k ) — ^ z —|- 4/i 3 // + 6/iz/ 2 ^ S§(X k ) 

— (p i + 6p 2 v + 2v 2 )S§(X k ) — (2 p 3 + 4pv)S7(X k ) — (2p 2 + v)S 3 (X k ) 


pS 9 (X k ) - -S 10 (X k ), 


(B 8 ) 
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Cl = -jU 4 - iw :i Si(X k ) - Q/xV + v 3 ^J S 2 (X k ) 

- (/j, 3 v + 3 fiu 2 )S 3 (X k ) - + 3/x 2 // + S 4 (X k ) - (/x 3 + 3/xi y)S s (X k ) 

- Q/x 2 + S 6 (X k ) - »S 7 (X k ) - |5 8 (X fe ), (B9) 


c 2 = -^ 3 -/xx/ 2 5'i(X fe )-(^ 2 ^ + xx 2 )S' 2 (X fc )-|(/x 2 + 6xx)5' 3 (X fe ) 

- (/x 2 + u)S 4 (X k ) - vS 5 (X k ) - ^S 6 (X k ). (BIO) 

The removal from a general quintic of the three terms in X 4 , X 3 and X 2 brings it to the 
Bring-Jerrard form (2.10), here re-written for convenience as 

X 5 + d x X + d 0 = 0. (Bll) 

By virtue of the Newton formulas (A3), the power sums for the quintic (Bll) are 


Si(Z k ) = S 2 (Z k ) = S 3 (Z k ) = 0, S 4 (Z k ) = —4di, S 5 (Z k ) = -5 d 0 . 


(B12) 


Assuming now, following Bring [49]], that the roots Z k of Eq. (Bll) are related by the quartic 
transformation (A8) to the roots Y k of the principal quintic (A7), we can substitute Eq. (A8) 
into Eq. (B12). This leads to a system of five equations with six unknown variables. More 
precisely, from the equation 


Si(Z k ) = 5xx 4 — 4ci — 3 xxiC 2 = 0, 


(B13) 


one finds 


The second equation [7 


4 3 

U 4 — -Cl + jC 2 U\. 


4 4g 

S 2 (Z k ) = -10XX1XX 2 C 0 - 4(lX 2 ) 2 Ci + -(Cl) 2 + 8C 0 C 2 + yXXiCiC 2 


6 


(xxi) 2 + 6 xx 2 


(c 2 ) 2 - 2 xx 3 (5c 0 + 4xxiCi + 3 u 2 c 2 ) = 0, 


obtained from the identities 


(B14) 


(B15) 


S 2 {Z k ) — Ss(Y k ) + 2uiSj{Y k ) + [(«i) 2 + 2u 2 \S§(y k ) + 2(u 4 u 2 + u 3 )S 3 (Y k ) 

+ [(^2) 2 + 2xx 4 + 2xx 1 xx 3 ]S , 4(Efc) + 2 (xx 2 xx 3 + xx 1 xx4)S' 3 (bfc) + 5 (xx 4 ) 2 , (B16) 
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S 6 (Y k ) = 3(c 2 ) 2 , S 7 (Y k ) = 7 Cl c 2 , S 8 (E fc ) = 8c 0 c 2 + 4( Cl ) 2 , (B17) 

relates m 2 and M 3 . The clever idea of the Bring-Jerrard method lies in choosing m 2 in such a 
way that the coefficient of M3 in Eq. (B15) vanishes. By inspection one finds immediately 


m 2 — 


5 Co 
3 c 2 


4 ci 


3 c 2 


Mi. 


(B18) 


Thus, Eq. (B15) now depends only on u\ and is a quadratic, i.e. |72] 


27(c 2 ) 4 - 160(d) 3 + 300c 0 dc 2 (m!) 2 + 27 Ci (c 2 ) 3 - 400c 0 (d) 2 + 375(c 0 ) 2 c 2 
+ 18(cic 2 ) 2 - 45c 0 (c 2 ) 3 - 250(c 0 ) 2 d = 0. 


Mi 

(B19) 


Last, by setting the sum of the cubes of (A8) to zero by virtue of (B12), a cubic equation 
for u 3 is obtained, by virtue of the identity 

12 

S 3 (Z k ) = 5 (u 4 ) 3 + ^ 6zSi(n), (B20) 

1 =2 

where (recall that we already know Si(Y k )...S 8 (Y k )) 


b 2 = 3m 2 (m 4 ) 2 , 

(B21) 

h = (M3) 3 + 3mi(m 4 ) 2 + 6m 2 m 3 m 4 , 

(B22) 

fr 4 = 3 (m 2 ) 2 m 4 + 3(m 4 ) 2 + 3m 2 (m 3 ) 2 + 6m 4 m 3 m 4 , 

(B23) 

h = 3(m 2 ) 2 m 3 + 3mi(m 3 ) 2 + 6m 4 (m 3 + m 4 m 2 ), 

(B24) 

b(j = (u 2 ) 3 + 3 (mi) 2 m 4 + 3(m 3 ) 2 + 6m 2 (m 4 + M1M3), 

(B25) 

b 7 = 3(mi) 2 m 3 + 3(m 2 ) 2 mi + 6(m 4 m 4 + m 2 m 3 ), 

(B26) 

b 8 = 3m 4 + 3(m 4 ) 2 m 2 + 3(m 2 ) 2 + 6miM 3 , 

(B27) 

b 9 = (mi) 3 + 3m 3 + 6MiM 2 , 

(B28) 

610 = 3m 2 + 3(mi) 2 , 

(B29) 

bn = 3mi, 

(B30) 

bi 2 = 1, 

(B31) 

Sg { Y k ) = 9coCi — 3 (c 2 ) 3 , 

(B32) 

Sio(bfc) = 5(c 0 ) 2 - 10 ci(c 2 ) 2 , 

(B33) 
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Su(Y k ) = -llc 0 (c 2 ) 2 - ll( Cl ) 2 c 2 , (B34) 

S V2 (Y k ) = -24 c oCl c 2 - 4( Cl ) 3 + 3(c 2 ) 4 . (B35) 


All intermediate quantities for reduction to the Bring-Jerrard form can be therefore found 
in terms of radicals. Of course, once the irreducible quintic (Bll) is solved in terms of hy¬ 
pergeometric functions as outlined in Sec. II, one has to invert the Tschirnaus-Bring quartic 
transformation (A8) to obtain the solutions of the principal quintic (A7) and, eventually, of 
the original quintic (Al), i.e., 

X k = X kiYjiZt)). (B36) 


Thus, one obtains in general twenty candidates for five solutions, and only numerical testing 
can tell which ones are correct 72]. 
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